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ABSTRACT 

We present the results of a search for cosmic ray point sources at energies above 4.0 x 10"'^^ eV in 
the combined data sets recorded by the AGASA and HiRes stereo experiments. The analysis is based 
on a maximum likelihood ratio test using the probability density function for each event rather than 
requiring an a priori choice of a fixed angular bin size. No statistically significant clustering of events 
consistent with a point source is found. 

Subject headings: cosmic rays — acceleration of particles — large-scale structure of universe 



1. INTRODUCTION 

The world data set of ultra-high energy cosmic rays is 
currently dominated by data recorded with the Akeno 
Giant Air Shower Array (AGASA) between 1989 and 
2003. AGASA has published arrival directions of 57 
cosmic-ray events above 4.0 x 10^^ eV obtained up to 
May 2000 ( Havashida et al. 2000). 

The High Resolution Fly's Eye (HiRes) experiment in 
Utah has operated in stereo mode since 1999. The pub- 
lished AGASA data set and the HiRes stereo data set 
overlap only a few months in time, but both detectors ob- 
serve approximately the same part of the northern sky. 
With 27 events above 4.0 x 10^^ eV recorded through 
January 2004, the HiRes data set now contributes sig- 
nificantly to the world ultra-high energy data set: thus 
it becomes increasingly interesting to search for possible 
point sources in this combined data set. 

In combining the two data sets, we need to account for 
the different errors on the individual cosmic-ray arrival 
directions as well as for the different background expec- 
tations in the two experiments. In addition, the analysis 
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method should be "unbinned," at least in the sense that 
any error from binning is much smaller than other errors 
in the data. This means that we do not define a fixed 
maximum angular separation for events to form a mul- 
tiplet. Rather, we search over the whole HiRes/AGASA 
sky for possible point source positions. 

For analyzing a data set combining events with very 
different errors, the maximum likelihood method is par- 
ticularly well suited. We perform a likelihood ratio test 
of the hypothesis that several events in the skymap come 
from a common source. In other words, we test whether 
any given position on the sky harbors a source which 
contributes > 1 source events to the data set. The 
likelihood of this hypothesis is compared to the null hy- 
pothesis Us — Q and this likelihood ratio is maximized 
using Us as a free parameter. By calculating the likeli- 
hood ratio for a dense grid of points on the sky, we essen- 
tially search the entire sky for the most likely position of 
a source of events. The statistical significance can be 
estimated by applying the same method to a large set of 
random isotropic data sets and evaluating what fraction 
of them have a likelihood ratio which is equal to or larger 
than the ratio observed in the real data described above. 

It should be noted that it is not obvious that a 
study of cosmic ray arrival directions will help iden- 
tify their sources, because charged particles suffer de- 
fiections in Galactic and extragalactic magnetic fields. 
The strengths of these fields are poorly known. Several 
authors have recently published estimates on the size of 
the expected deflections, and the results differ c onsid- 
erably ((Sid. Miniati & Ensslln 2004; Dolag et al.ir2003) . 
The point source search reported in this paper could be 
fruitful only if deflections at energies above 4.0 x 10^^ 
eV are sufficiently weak, making cosmic ray astronomy 
possible. 

The AGASA data set has b een intensively studied 
and described in t he literature llHavashida et al. 1 119961: 
iTakeda et al.F^IHavashida et al.l200nHTa,keda et al.l 
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|2003)- The HiRes stereo data set is relatively new, so 
we will describe it in more detail in Section 2. Section 3 
gives a description of the likelihood method applied here. 
Results from the analysis of the combined data set are 
described in Section 4. In Section 5, we test the likeli- 
hood method with simulated data sets, and a discussion 
follows in Section 6. 

2. THE HIRES DETECTOR 

HiRes is a stereo air fluorescence experiment with two 
sites (HiRes 1 and 2) at the US Army Dugway Prov- 
ing Ground in the Utah desert (112° west longitude, 
40° north latitude, with a vertical atmospheric depth of 
860g/cm^). The two sites are separated by a distance of 
12.6km. 

The ultra-high energy cosmic ray flux is small and is a 
steeply falling power law in energy. Thus experiments at 
ultra-high energies need a large detector volume. Con- 
sequently, the primary cosmic ray particles can not be 
observed directly, since they interact in the upper atmo- 
sphere and induce extensive air showers with of the order 
of 10^" particles for a 10^^ eV primary. The properties 
of the original cosmic ray particle, such as arrival direc- 
tion and energy, have to be inferred from the observed 
properties of the extensive air shower. In HiRes, this 
is achieved by observing the fluorescence light produced 
when particles of the extensive air shower interact with 
nitrogen molecules in the atmosphere. This method has 
the advantage that the shower development in the at- 
mosphere is imaged and important quantities like the 
shower size and the height of the shower maximum can 
be determined directly. If the shower is viewed simul- 
taneously by two detectors in stereo mode, the arrival 
direction can be reconstructed with an accuracy of less 
than 1°. The main shortcoming of the technique is the 
low duty cycle of only about 10%, as air fluorescence 
detectors can only be operated on dark, moonless nights 
with good atmospheric conditions. 

To observe fluorescence light from air showers, the de- 
tector at each site is made up of several telescope units 
monitoring different parts of the night sky. With 22 (42) 
telescopes at the first (second) site, the full detector cov- 
ers about 360° (336°) in azimuth and 3° -16.5° (3° -30°) 
in elevation above horizon. Each telescope consists of a 
mirror with an area of about 5 m^ area for light collec- 
tion and a cluster of 256 photomultiplier tubes in the 
focal plane. 

The shower geometry is determined by a global min- 
imization using both the timing and pointing informa- 
tion from all tubes. From measurements of laser tracks 
and stars in the field of view of the cameras we estimate 
that the systematic error in the arrival direction deter- 
mination is not larger than 0.2°, mainly caused by un- 
certainties in the survey of mirror pointing directions. 
Various aspects of the HiRes detector and the recon - 
struction procedures a r e described in Bove r et al.l l)2002|) : 
iSadowski et"aD l|2002|) : lMatthews et al.l (2003). 

For the present analysis, we use stereo data taken be- 
tween December 1999 and January 2004. This sample 
is subje ct to the same quality cuts used in I Abbasi et al.l 
l)2004bt) . A minimum track length of 3° in each detector, 
an estimated angular uncertainty in both azimuth and 
zenith angle of less than 2° , and a zenith angle less than 
70° are required. In addition, the estimated energy un- 



certainty is required to be less than 20% and x^/dof < 5 
for both the energy and the geometry fit. A total of 27 
events above 4.0 x 10^^ eV pass the selection criteria. 

An atmospheric data base built from the reconstruc- 
tion of laser shots is used to correct for the impact of the 
atmospheric conditions on the reconstruction of events. 
The data base contains hourly values of parameters de- 
scribing the aerosol content of the atmosphere. Simula- 
tions show that after correcting for the aerosol content, 
the impact on the event geometry and energy is small for 
a wide range of different atmospheric conditions. 

The angular resolution of HiRes is determined 
using simulated showers generated with CORSIKA 
6 (Heck ct al. 1998) and QGSJET for the first interac- 
tion. After the quality cuts described above, 68 % of all 
showers generated at 10^^ eV are reconstructed within 
0.57° of the true shower direction. The angular reso- 
lution depends weakly on energy. The 68 % error ra- 
dius grows to 0.61° and 0.69° for showers generated at 
4.0 • 10^^ eV and 10^*^ eV, respectively, because showers 
at higher energy are farther away, on average. 

3. ANALYSIS 

The maximum likelihoo d method used here is outlined 
in lKinnison et al.l ljl982() : a useful and clarifying appli- 
cation to the problem of finding a deficit from the posi- 
ti on of the moon in a cosmic ray skymap may be found 
in Wascko' (2000") . For a general description of likelihood 
methods, see for example Eidelman et al.. (.2004) : Meverl 
( 1975) and references therein. 

The likelihood method can be applied both in searches 
for emission from a priori selected source locations, and 
in searches for the most likely location of such a point 
source on the whole sky. 

Consider first a fixed source location in right ascension 
and declination: Xg = (as,^s)- Given a sample of N 
cosmic ray events, we suppose that Ug events come from 
the source location Xg, and that the remaining N — rig 
events are random background events. If the ith event 
is a source event, then its true arrival direction is Xg. 
The probability density for finding it at some location x 
is given by the function Qi(x,Xg), derived from the an- 
gular errors of the event and the angular displacement 
between Xg and x. On the other hand, if the ith event 
is a background event, then the probability density for 
finding it at a location x is given by the function Ri(x), 
derived from the relative exposure of the detector to an 
isotropic background of cosmic rays. The subscript is 
necessary because R depends on whether event i is a 
HiRes or AG AS A event. Each of these functions is nor- 
malized to unity over positions x in the sky. 

We do not hypothesize which individual events are 
source or background events. We only suppose that there 
are rig events in the sample that come from some source 
position Xg. Therefore, the partial probability distribu- 
tion of arrival directions x for the ith event is given by: 

N ~ Ug 

R,{x,xg) ^ —Qi{x,Xg) ^ — — Ri[x) . (1) 

It follows that the probability of finding the ith event 
at the location Xi (where it is actually observed) is 

Pi (^i 5 ) ■ 

The likelihood for the entire set of N events is then 
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given by: 



N 



(2) 



i=l 



The best estimate for the number of source events, un- 
der the assumption of a point source located at xl, is 
determined by finding the value of rig which maximizes 
C. 

Because the value of the likelihood function depends 
on the number of events, a more useful quantity than C 
is the likelihood ratio TZ: 



C{0,xs) 



(3) 



where C{0,Xs) is the likelihood function of the null hy- 
pothesis (ris = 0). In practice, we maximize InTZ, which 
is equivalent to maximizing C. 

The method described so far is sufficient for testing a 
specific source position Xs- To search the entire sky for 
the source position with the strongest signal, we calculate 
the likelihood ratio hiTZ for a dense grid of points on 
the sky covering the full range of equatorial coordinates 
accessible to AGASA and HiRes. The source position 
is essentially treated as a free parameter, along with the 
number of source events Ug . Searching for the parameters 
as, Ss, and which maximize the likelihood ratio will 
therefore give us the best estimate for the position of the 
source and the number of events it contributes. 

The search proceeds as follows. The visible region of 
the sky is divided into a fine grid of points with separa- 
tions 0.1°/ cos (5 and 0.1° in a and S, respectively. Each 
point is treated in turn as a source location Xs- At each 
point, the specific quantities Qi(xi,Xs) and Ri{xi) are 
required for every event. The source probability density 
function Qi{xi, Xg) depends on the angular resolution as- 
sociated with the zth event, which in principle may de- 
pend on several quantities including energy, zenith angle, 
etc. The background probability density function Ri{xi) 
depends on the detector exposure to different parts of 
the sky: it is generally the same function for all events 
observed by a given detector, but may in principle be a 
function of e.g. energy as well. For each event, Qi{xi, Xg) 
needs to be reevaluated at every grid point, while Ri{xi) 
needs only to be evaluated one time. 

Once the values of Qi and Rt are specified, the log 
likelihood ratio IxiTZ. is maximized with respect to ris, 
where > 0. This process is repeated for each position 
Xs on the grid. For most locations. In 7^ is zero, but local 
maxima will occur in the vicinity of one or more events. 
Note that while the method is in some sense binned due 
to the discrete array of grid points, the spacing is chosen 
small enough so that errors introduced by binning are 
insignificant. 

4. RESULTS 

The procedure described above is applied to the com- 
bined data set of 57 AGASA events and 27 HiRes stereo 
events above 4.0 x 10^^ eV. The analysis is restricted to 
the AGASA field of view with -10° < ,5 < 80°, and aU 
but one HiRes event fall into this declination range. 
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Fig. 1. — Normalized probability densities (p.d.f.) for random 
background in HiRes and and AGASA, in equatorial coordinates. 
The sky is binned in 0.1°/ cos 5 and 0.1° bins in right ascension a 
and declination 5, respectively. 



The background expectations for HiRes and AGASA 
are different. The HiRes exposure has some dependence 
on right ascension. Fig.^ shows the normalized HiRes 
and AGASA background map in equatorial coordinates. 

For the signal probability density function (Q) of the 
HiRes events, we use a two-dimensional Gaussian func- 
tion whose width is chosen such that 68 % of the proba- 
bility density function falls within an opening angle 0.6°. 
Note that for a two-dimensional Gaussian distribution 
the opening angle = 1.515 ti encloses 68 % of the distri- 
bution. Since the dependence on energy is weak, we use 
the same value, a = 0.4°, for every HiRes stereo event. 

For AGASA, we approximate the probability density 
by the sum of two Gaussian functions ch osen such that 
the 68 % and 90 % opening angle given in iTakeda et al.l 
l)1999(l is correctly reproduced. 
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The width of the Gaussians as a funct ion of energy 
has the same energy dependence used in lAbbasi et al.l 
lf20f)4^ . and is given by 

ai = (6.52° - 2.16° log(£;Ecv))/l-515 and (5) 
(72 = (3.25°- 1.22° log(^Ecv))/1.515 . (6) 

More detailed information on the error shape of indi- 
vidual AGASA events could easily be implemented if it 
becomes available. 

Fig-El shows the result of the analysis. At each a and 
5, the likelihood ratio is shown for the number of source 
events Us which maximizes In TZ. One can clearly recog- 
nize where events are located, and one can also recognize 
locations with several nearby events. AGASA and HiRes 
events can easily be distinguished, as the latter have bet- 
ter resolution and therefore smaller regions with large 
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Fig. 2. — Likelihood ratio InTi, maximized with respect to Us, as a function of right ascension and declination of the source position for 
the combined set of AGASA and HiRes events above 4.0 X 10^^ eV. Local maxima occur wherever events or clusters of events are located on 
the sky. The global maximum, i.e. the most likely position of a "point source" is at right ascension a = 169.3° and declination 5 = 57.0°. 
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Fig. 3. — Logarithm of the likelihood ratio as a function of the 
number of source events for the position of the maximum. 
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Fig. 4. — Maximum InTZ for 10* simulated random data sets 
with the same number of AGASA/HiRes events as the actual data 
set. The hatched area marks random sets whose maximum InTZ 
exceeds the value for the real data set. 



likelihood. The point with the largest InTZ is at right 
ascension a — 169.3° and declination S = 57.0°. The 
corresponding event cluster comprises 3 nearby AGASA 
events with coordinates {a, S) and energies E of (1) 
(168.5°, 57.6°), E = 77.6 EeV, (2) (172.3°, 57.1°), E = 
55.0 EeV, and (3) (168.3°, 56. 0°), E = 53.5 EeV. This 
cluster has been describe d in iTakeda et al.l lll999D and 
is listed as cluster C2 in iHavashid^^T^Ll l|2000() . The 
maximum likelihood ratio at this position is In??, = 8.54 
for Hs = 2.9. 

The behavior of In 7^ as a function of the fit parameters 
(number of source events and position of the source) in 
the vicinity of the maximum gives error estimates for the 
fit parameters. To a good approximation, the 1 a error is 
given by the interval over which InTZ drops by 0.5 from 
its maximum value. 

Fig.O shows In??, as a function of the number of source 
events at the position of the maximum. The best es- 
timate of the number of source events is ris = 2.9^^1- 
Similarly we find a = 169.3° ± 1.0° and 6 = 57.0° ± 0.5° 
as the best estimate for the position of the maximum. 

The statistical significance of the appearance of a 



"source" with a maximum likelihood ratio InTZ in the 
combined AGASA/HiRes data set can be evaluated us- 
ing simulated random data sets. The full likelihood anal- 
ysis is performed for 10* random data sets with the same 
number of AGASA/HiRes events and the same underly- 
ing exposure as the original data set, but isotropic ar- 
rival directions. The chance probability for the "source" 
to appear is then given by the fraction of random data 
sets which have at least one location causing the maxi- 
mum In 7^ to be equal or larger than 8.54, the value of 
the maximum in the real data. 

Fig.0] shows the distribution of the maximum InTZ for 
each these random data sets. Out of 10* simulated data 
sets, 2793 have a maximum InTZ exceeding that of the 
real data set. The chance probability of the source hy- 
pothesis is therefore of the order of 28%. Consequently, 
there is no statistically significant evidence for clustering 
consistent with a point source in the combined data set. 

Note that this is not simply the chance probability for 
a triplet, but rather the chance probability for a set of 
27 HiRes events and 57 AGASA events to contain a "hot 
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Fig. 5. — Top: Maximum In 7^ for 10* simulated random data 
sets where a 4-event source (blue line) or a 3-event source (red line) 
has been added. Bottom: Fraction of simulated data sets with 
maximum InR. > InlZcut as a function of InTZcut for simulated 
random isotropic data sets with no source added (black line), a 
4-event source added (blue line), and a 3-event source added (red 
line). 



spot" with as high a probabihty to be a "source" as the 
triplet. Many of the simulated likelihood ratios larger 
than 8.54 in Fig.^jare indeed caused by doublets. 

5. TEST OF THE METHOD WITH SIMULATED DATA SETS 

The maximum likelihood method is tested by applying 
it to simulated data sets with sources. These simulated 
data sets have m events from a common source added 
to an otherwise isotropic arrival direction distribution. 
To create such a source, we pick a point in the sky for 
the source location and generate m events with arrival 
directions deviating from the source location according 
to the probability density function of the individual event 
as described in Section 2. The source location is chosen 
at random, but the distribution of locations is forced 
to reflect the overall exposure of the detector (Fig.P), 
so that regions with higher exposure are more likely to 
contain a source. The m source events replace events 
in the data set, so the total number of events is always 
83. There might be additional events close to the source 
location due to chance. 

The full likelihood analysis is applied to these random 
data sets. Fig.El(top) shows the maximum InTZ for 10^ 
random data sets where a source with m — A (blue line) 
or m = 3 (red line) has been added. Fig.|Sl(bottom) 
shows the fraction of simulated data sets with maxi- 
mum InT^ > InTZcut as a function of InTZcut for ran- 
dom isotropic data sets with no source added (m = 0, 
black line), a 4-event source added (blue line), and a 
3-event source added (red line). While there is sub- 
stantial overlap between the distributions, the plots also 
show that point sources which add three or more events 
to an isotropic map are recognized with high efficiency. 
The medians of the three maximum In 7?.-distributions 
(to = 0,3,4) are 7.7, 11.6 and 16.0, respectively. 

The HiRes and AGASA events in the sample have 



number of HiRes events in 3 (4) source events 

Fig. 6. — Median of the maximum In 7^-distribution for simulated 
data sets with m = 3 and m = 4 source events as a function of the 
number of HiRes events contributing to the source cluster. 
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Fig. 7. — Difference between the fitted and true right ascension 
a and declination S for 10** data sets with an artificial point source 
added to an isotropic data set for m = 3 (left) and m = 4 (right) 
source events. The width of a Gaussian fitted to the distribution 
is indicated. 



rather different angular errors, and we expect this to be 
reflected in the likelihood ratio of the simulated data sets. 
The maximum In 7?.-distribution for clusters that are 
dominated by HiRes events should have a larger median 
than the data sets dominated by AGASA events. Fig.El 
shows the median of the maximum In 7?.-distribution as 
a function of the number of HiRes events that contribute 
to the cluster. Both for data sets with to = 3 and to = 4, 
the median increases if more HiRes events are part of the 
cluster. 

Artificial sources can also be used to test the accu- 
racy with which the position of the point source is recon- 
structed. Fig.[71 shows the difference between the fitted 
and the true right ascension a and declination angle S for 
TO = 3 and to = 4. The distributions are fit to a Gaus- 
sian function, and the width, averaged over all source 
locations and combinations of AGASA and HiRes source 
events, is of order <t = 0.3° to cr = 0.5°. The resolution 
improves with the number of events contributing to the 
signal, as expected. 

6. DISCUSSION 

The chance probability of th e triplet using AGA SA 
data alone has been estimated in lTakeda et al.l l)2003(l as 
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being of order 1 %. This estimate is based on a fixed bin 
size of 2.5°. To test what chance probabihty an unbinned 
analysis gives, we repeat the likehhood analysis for the 
data set comprising only the 57 AGASA events above 
4.0 X lO^^eV. The largest likelihood ratio {InTZ = 9.66) 
appears again near the events forming the triplet, with 
a = 169.3° ± 1.0° and S = 57.0° ± 0.5° for n, = 2.9tl°. 

As before, we evaluate the chance probability for the 
appearance of a source with maximum InTZ — 9.66 or 
higher in this data set by analyzing a large number of 
simulated isotropic data set, now containing 57 AGASA 
events. 452 out of lO'' random data sets have a maxi- 
mum In 72. in excess of 9.66, so the chance probability is 
4. 5%. The increase ov er the chance probability given 
in iTakeda et al. I (|2003^ reflects the fact that the un- 
binned maximum likelihood analysis evaluates the qual- 
ity of the cluster on a continuum. This improves the 
sensitivity to a true point source, and at the same time 
it removes the artificial fluctuations that arise in a binned 
analysis when an event falls just inside or outside the bin. 

The maximum likelihood method allows a straightfor- 
ward evaluation of the chance probability that a given 
cluster of events come from a common source rather 
than being caused by random background. It can also 
be used to estimate the chance probability that one or 
more pre-defined positions on the sky correlate with 
cosmic ray arrival directions. If several sources con- 
tribute, but none of them are strong enough to pro- 
duce a notable value of In TZ, then the likelihood method 
described here has to be complemented by other meth- 
ods, such as the angular two-point correlation function. 
Searches for excess clustering have been performed us- 
ing both the AGASA and HiRes data sets separately, 
and for the combined set. N o significant clusteri ng is 
reported in the HiRes stereo (* Abbasi et al.ll2004bD and 
monocular data sets (Abbasi et al. 2004c), and earlier re- 
port s that the AGASA d ata set shows significant cluster- 
ing l)Takeda et al. llT999tl have re centlv bee n ques tioned. 
An unbiased analvsis ijFinlev fcW cstcrhoffl l2004| ) which 
sets aside the fraction of the AGASA data that was 
used to formulate the clustering hypothesis, finds that 
the clustering is consistent at the 8 % level with the null 
hypothesis of isotropically distributed arrival directions. 
In addition, a two-point correlation analysis of the com- 
bined AGASA and HiRes data sets finds no significant 
evidence for clustering (Finlcv et al. 2004). 

The energy threshold of 4.0 x 10^^ eV chosen in this 
analysis is dictated by the fact that only AGASA data 
above this threshold is published. I t has been sug- 
gested (e.g. iDe Mar co. Blasi. fc Ohnto' (20031)) that the 
discrepancy in the m easurement of the cosmi c ray flux 
by AGASA ( T akeda et al. 1998; Havashida et al. 200(| 

and the monocular HiRes detector (Abbasi et al 2004afl 

can be reconciled to some degree by an ad hoc change 
of about 30 % in the relative energy scales of the ex- 
periments. In this scenario, numerically equal energy 
thresholds in HiRes and AGASA are not equivalent, and 
in order to compare the two data sets, one should either 
raise the nominal AGASA threshold or lower the HiRes 
energy threshold, or apply a combination of both shifts. 

We do not speculate on the nature of the flux discrep- 
ancy here. At this point, there is no indication that it is 
caused by such a systematic shift in the energy scale of 
one or both of the experiments. We also note that lower- 



ing the energy threshold could potentially make magnetic 
bending a more serious limitation. However, since the 
possibility of an energy shift has been suggested by oth- 
ers, it is useful to evaluate whether the conclusions of the 
maximum likelihood analysis change under this assump- 
tion. Therefore, the method is applied again to the data 
set combining the 57 AGASA events above 4.0 x 10^^ eV 
and the 40 HiRes events above 3.0 x 10^^ eV, an energy 
threshold roughly 30 % lower. 

Lowering the HiRes energy cut admits an event that 
is near the aforementioned AGASA triplet. This HiRes 
event has an energy of 37.6 EeV and is located at right 
ascension a — 169.0° and declination 5 = 55.9°. The 
location with the largest \nTZ is close to this event, at 
a = 169.1° ± 0.6° and 6 = 56.3° ± 0.4°. The maximum 
likelihood ratio is In 72 = 12.98 for = i-Qt^i^- 47 out 
of 10* random isotropic data sets have a larger value of 
maximum In TZ. 

It is not possible to convert this number into a valid 
chance probability for the combination of the AGASA 
triplet and the additional HiRes event because it contains 
several biases. While the maximum likelihood method 
removes the bias due to choosing an angular bin size, 
it does not remove the bias due to choosing an energy 
threshold. The AGASA th reshold of 4.0 x lO^^ eV was 
chosen in iHavashida et al.l (1996) precisely because it 
maximized the clustering signal in the early data set, 
a data set which comprises the first 30 of the 57 AGASA 
events (and two of the triplet events). An unbiased 
search for point sources using this energy threshold would 
require that these 30 events be set aside. Further bias 
occurs because of the ad hoc lowering of the HiRes energy 
threshold. 

The best estimate for the position can be considered 
an a priori location for a point source search with sta- 
tistically independent future data. To illustrate the sen- 
sitivity of a new data set to this hypothesis, simulations 
show that a simple binned analysis with a 1° source bin 
radius centered on this location would give an unbiased 
chance probability of 10~^ if 2 out of 40 future events 
fall into the source bin. In this calculation, it is assumed 
that the geometrical acceptance for the 40 new events is 
the same as for the original 40 events. 

HiRes is currently the only detector observing the cos- 
mic ray sky in this energy range from the northern hemi- 
sphere. We anticipate several more years of data taking, 
and the search for cosmic ray point sources at the highest 
energies will continue when more data becomes available. 
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